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Fast magnetic reconnection refers to annihilation or topological rearrangement of magnetic fields on 
a timescale that is independent (or nearly independent) of the plasma resistivity. The resistivity of 
astrophysical plasmas is so low that reconnection is of little practical interest unless it is fast. Yet, 
the theory of fast magnetic reconnection is on uncertain ground, as models must avoid the tendency 
of magnetic fields to pile up at the reconnection layer, slowing down the flow. In this paper it is 
shown that these problems can be avoided to some extent if the flow is three dimensional. On the 
other hand, it is shown that in the limited but important case of incompressible stagnation point 
flows, every flow will amplify most magnetic fields. Although examples of fast magnetic 
reconnection abound, a weak, disordered magnetic field embedded in stagnation point flow will in 




general be amplified, and should eventually modify the flow. These results support recent arguments 
against the operation of turbulent resistivity in highly conducting fluids. © 1998 American 
Institute of Physics. [S1070-664X(98)04801-0] 


I. INTRODUCTION 

Magnetic fields in astrophysical systems are almost com- 
pletely frozen to the plasma, because the magnetic Reynolds 
number R m , the ratio of the Ohmic decay time to the dy- 
namical time, is extremely large; of order 10 15 to 10 2! for 
interstellar fields and 10 8 to 10 10 for stellar fields, for ex- 
ample. Any departure from frozen in behavior — i.e. any re- 
connection of the magnetic fieldlines — is a finite conductiv- 
ity effect. If magnetic reconnection occurs on a timescale 
which is independent of R m , the reconnection is said to be 
fast. Petschek’s model 1 at its maximum rate is almost fast, 
depending only logarithmically on R m . The Sweet and 
Parker time independent model 2,3 is slow, with the reconnec- 
tion rate scaling as R m 1/2 . Linear tearing modes 4 grow at a 
rate scaled by R~ V5 and are also slow. 

Given the large value of R m in many astrophysical situ- 
ations, reconnection can be of little practical importance un- 
less it is fast. Therefore, the problem of developing models 
for fast magnetic reconnection is of considerable interest. 

Clark 5 and later Moffatt 6 developed an analytical model 
for almost fast magnetic reconnection based on two dimen- 
sional (2D), incompressible flow near a hyperbolic stagna- 
tion point. I refer to this model as 2DHS throughout the text. 
The magnetic field is initially amplified by induction, and 
then decays resistively at a superexponential rate. Strauss 7 
showed that this simple model describes a numerical simu- 
lation of magnetic reconnection in a current sheet. 

The 2DHS model is kinematic, in the sense that the flow 
field is prescribed. During the transient amplification phase, 
magnetic forces can become very large, scaling as and 
one concludes that the reconnection is quenched or greatly 
modified by these large forces, even though they act for only 
a short time. The dynamical, steady state reconnection model 
of Craig and Henton 8 shows magnetic flux pileup and large 
pressure gradients near the stagnation point of the flow — 
which is also an X-point of the magnetic field — and might 
represent the outcome of the 2DHS model if magnetic forces 
were accounted for. See References 5, 9 and 10 for other 


discussions of the role of large pressure gradients in sustain- 
ing the flow in fieldline reconnection. 

The 2DHS model exemplifies a generic difficulty of re- 
connection theory: in order that resistive effects occur rap- 
idly, the magnetic field must develop structure on small spa- 
tial scales which are proportional to a positive power of the 
plasma resistivity 77. These small spatial scales generally im- 
ply large Lorentz forces, scaling as a negative power of rj. 
But, fast reconnection requires that the plasma velocity re- 
main independent of tj. 

The difficulties associated with the 2DHS model can be 
made to disappear in three dimensional flows. More gener- 
ally, it turns out to be relatively easy to characterize the 
action of any linear stagnation point flow on any magnetic 
field, and to set forth conditions under which the fields decay 
and Lorentz forces do not grow, so that if the kinematic 
reconnection theory is self-consistent initially then it remains 
so. This allows one to estimate, for example, the degree of 
reconnection of a weak magnetic field in a turbulent fluid. It 
may also be a useful preliminary step toward developing 
fully dynamical models of reconnection, which apply when 
the fields are strong. 

In Sec. II of this paper I review the 2DHS model and 
generalize it to three dimensions. In Sec. Ill I develop the 
solution for arbitrary fields in arbitrary hyperbolic stagnation 
point flows. In Sec. IV D I summarize the implications of 
this work for astrophysical reconnection, and for the opera- 
tion of turbulent resistivity, which is often invoked in astro- 
physics. 


II. THE 2D HYPERBOLIC STAGNATION POINT 
MODEL AND EXTENSION TO THREE DIMENSIONS 

Consider a magnetic field B= (#(y, 0,0,0) embedded in 
a bulk flow u with a hyperbolic stagnation point at the origin, 

u=w'(*,-y,0). (I) 
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Equation (1) should of course be regarded as a local repre- 
sentation of a flow which is globally bounded. Assume the 
medium has resistivity rj and Ohmic diffusivity \ 
= t}c 2 /4tt. The magnetic induction equation, 

SB 

— = Vx(vxB) + \V^B, (2) 


reduces in this case to 

SB SB S 2 B 

— =u'v — b U ' B + A. — y , 

St * Sy Sy 

Taking as the initial condition 

B(y,0) = B 0 sin k 0 y. 


(3) 

(4) 


Eq. (3) has the solution 

B(y,t) = B 0 e u ''- R » {e2u '-' ) sin k 0 e u ''y , (5) 

where the magnetic Reynolds number R m =(\kl/2u ' ) is 
roughly the ratio of the initial Ohmic decay time to the flow 
time within one magnetic spatial period. (Although in steady 
state problems it is standard procedure to define the dynami- 
cal time from the advective term in the induction equation, 
so that R m is linear in the magnetic lengthscale, the definition 
used here is natural in the present problem.) It is the depen- 
dence on or 77, which is most important here. Equation (5) 
describes a field which is initially amplified by compression 
as it is swept toward the x axis, while the scale over which 
the field reverses shrinks exponentially. Eventually Ohmic 
processes dominate compression and the field begins a phase 
of superexponential decay. The fieldstrength reaches its 
maximum, B max , at time t max , 


t max 


2 u’ 


1 

Ini -R„ 


B n 


Bn 



1/2 

sin k 0 1 



(6) 


where the approximation consists of assuming R m > 1. Thus, 
the field is amplified by a factor of order R^ 2 or 77 ~ l/2 , and 
the magnetic lengthscale decreases by the same factor. The 
decay phase is so fast that the field falls to Me of its initial 
value in a time of order (2 « f ) _l ln(/? m In R^ 2 ). In view of 
the weak dependence of the decay time on R m , this system is 
a model of “almost fast” magnetic reconnection. However, 
it is easy to show that the current density and Lorentz force 
also peak before decaying, and that the time integrated Lor- 
entz force F 1 is of order 


/ 


dtF L ~ 


4 7TM ' 


(7) 


The fluid will be slowed down if the time integrated Lorentz 
force is greater than or equal to the initial momentum density 
in the fluid, measured, say, at a distance k 0 1 from the origin. 
Therefore, the fluid is decelerated unless B 0 
< v(4 7rp)(w '/&())/?„, 3/4 , which means that the initial Alfven 
Mach number of the flow must exceed fi 2/4 . Thus, this re- 
connection model is self-consistent only for extremely weak 
magnetic fields. It is difficult to imagine applying it in un- 


modified form to the interstellar medium, where B is close to 
equipartition with the turbulent gas velocity. Similar objec- 
tions would arise if the model were applied to the convection 
zones of stars; even more so in stellar coronae, which are 
magnetically dominated. 

These difficulties can be avoided in a three-dimensional 
flow'. Consider the velocity field, 

u=«'(— *,-y, 2 z). ( 8 ) 


Fluid is swept in along the x and y axes and ejected along the 
Z axis, and the flow is incompressible (see Reference 1 1 for 
a discussion of dynamically consistent flows of this type). 
Again, assume the magnetic field points in the a direction 
and depends only on y and t. The induction equation (2) is 

SB SB S 2 B 


— = w ' y - — u f B + \ 
St ' Sy 



(9) 


The solution which satisfies the initial condition (4) is 


B = B 0 e~ u ’' ^ ](e2u '~' ) sin V“V (10) 

Equation (10) differs from Eq. (5) in that the magnetic field 
is not compressed in the initial nonresistive phase; instead, 
the field is weakened by ejection in the z direction even as its 
spatial scale decreases. Once resistivity becomes important, 
which happens at about the time t max given in Eq. (jS), the 
spatial scale has decreased to R m 1/2 of its initial value, and 
decay is superexponential. The current density is nearly flat 
with time in the nonresistive phase, and the Lorentz force 
decreases monotonically, so if F L is initially negligible it 
remains so. 

The results of this section suggest that one could system- 
atically characterize magnetic reconnection near stagnation 
points in three dimensional flows, and select the flows and 
fields for which reconnection is fast and unaccompanied by 
large increases in the Lorentz force. This is taken up in Sec. 
III. 


III. GENERAL STAGNATION POINT FLOW 

The resistive induction equation cannot be solved ana- 
lytically for general flows, even when the velocity is a linear 
function of the coordinates. The solution to the nonresistive 
induction equation can, however, be written down exactly, 
and it is then easy to see whether the magnetic field develops 
progressively smaller scale structure over time, and how the 
Lorentz force evolves during the ideal phase. 

Much of this section depends on results from linear al- 
gebra, and are given in any one of a number of books. The 
author used the text by Curtis, 12 which includes material on 
differential equations. 

It is convenient to adopt a Lagrangian description of the 
flow. Label points in the fluid by their initial positions x 0 , 
and write their positions at subsequent times / >0 as x(x 0 ,/). 
One can then form the deformation matrix D at every point 
and time, which has elements 
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The volume element d*x changes with time according to 
d*x=\D\d*x 0 , where |Z)| is the determinant of D. The in- 
duction equation (2) with A = 0 can then be integrated to give 
B(x,0 in terms of the initial field B 0 , 


B(x,/) = 


D B Q (x Q (x,r)) 

\D\ 


( 12 ) 


Equation (12) shows that growth of the field is associated 
with an increase of the deformation matrix with time, which 
for incompressible or nearly incompressible flows is associ- 
ated with shearing or stretching of volume elements. The 
change in the magnetic lengthscale is also brought about by 
shear, as reflected in the map of x back to x 0 . In general, this 
map (or its inverse, the trajectory of a fluid particle) is diffi- 
cult to construct, but it becomes simple near a stagnation 
point, taken to be the origin of the coordinate system, at 
which the Lagrangian equations of motion are 


dx 

— = U x, 
dt 


(13) 


where U is a 3X3 matrix with constant, real coefficients. 
The solution of Eq. (13) with the initial condition x(0) = x 0 is 


x(t) = e' lu -x 0 , 

where the matrix e' v is defined by its Taylor series. 

04) 

.--tai ( f - 

n~+xj - 0 J • 

(15) 

It is clear from Eqs. (11) and (14) that 


D=^ rU , 

(16) 

and that the inverse of Eq. (14) is 


Xo(x,r) = e~' u -x. 

(17) 


Assume for the remainder of the paper that \e fU | = 1 (incom- 
pressible flow). Using this together with Eqs. (14) and (17), 
Eq. (12) becomes 

B(x,/) = e' u B 0 (<>-' lJ .x). (18) 

The matrix e' 11 is generally not easily computed as it stands. 
However, by a theorem of linear algebra, the matrix U is 
similar to the sum of a diagonal matrix M and a nilpotent 
matrix N which commutes with M, 

U = S(M+N)S -1 , (19) 

for some nonsingular matrix S. Because M and N commute, 
this carries over to exponentiation, 

e tV =Se t{M+ ^S~ l . (20) 

The elements of the diagonal matrix M are just the roots a ( 
of the minimal polynomial of U (the lowest degree polyno- 
mial equation satisfied by U), with each root occurring ac- 
cording to its multiplicity. It is readily shown that e m is 
diagonal, with element (e tM ) h = Computation of c' N is 
also straightforward; since N is nilpotent, the Taylor series 
given in Eq. (15) terminates at or before the N 3 term. 
Returning to Eq. (18), written now in the form 

B(x,f) = Se /(M + N) S _1 * B(Se _ ' <M+N) S~ 1 *x), (21) 


we can make the following observations. The matrix S is 
independent of time and does not directly cause growth or 
decay of the magnetic field or current. The factors e : " /N lead 
at most to algebraic growth or decay of the field and current, 
by virtue of the finite Taylor series which define these ma- 
trices. Exponential growth or decay of the field and its 
lengthscale comes about only through the action of the e~ m 
matrices, and the rates of growth/decay are easy to read off; 
they correspond to the real parts of the roots a, of the mini- 
mal polynomial of the original velocity matrix U, which are 
also eigenvalues of U. 

There is a constraint on the a, which arises from the 
incompressibility of the flow. It follows from Eq. (20) that 
|* ,M e fN |==|e ,M ||e ,N |:= 1. On the other hand, 

|e' M | = exp| <2 «,«,■). (22) 

where n is the multiplicity of the i th root. Since |e' N | is an 
algebraic function of f, the real part of the sum in Eq. (22) 
must be zero. 

In order for the lengthscale of the field to shrink to the 
resistive scale, at least one a, must have a negative real part. 
The field grows because of the action of the a t with positive 
real part. Because of the remark following Eq. (22), if there 
is shrinking there must also be stretching, so magnetic fields 
with arbitrary direction and coordinate dependence are am- 
plified even as their lengthscale shrinks. 

The velocity field given by Eq. (8) illustrates these con- 
clusions. In this case, the U matrix is diagonal. The magnetic 
field derivatives with respect to .r and y grow as e u The x 
and v components of B shrink by the same factor. The z 
component of B grows, however, as e 2u '. The Lorentz force 
is bounded with time only for fields in the x — y plane, so 
only fields confined to the x — y plane can undergo fast, ki- 
nematic reconnection. 

A simple modification of the flow (8) leads to algebraic 
growth of B. Consider the velocity field 

u=u’(-x-y,-y,2z). (23) 

In this case, direct integration of the equations of motion (13) 
yields 


x=x 0 e “ '-y 0 u'te " y=y 0 e z = Z 0 e 2u '', 


(24) 


while the matrix U is the direct sum of a diagonal matrix 
with elements (- 1,- 1,2) and the nilpotent matrix N with 
single nonzero element jV 12 =- 1, which satisfies the equa- 
tion N 2 = 0. In this example, the U matrix has a doubly de- 
generate eigenvalue - 1, and this corresponds to forcing of x 
by y at its natural decay rate. This is the origin of the com- 
bined exponential and algebraic deformation seen in Eq. 
(24). 

As an example of the action of this flow on a magnetic 
field, let 


B o(*o) = /(*o)(0.1.0)- (25) 

Using Eq. (12), the field at time I is seen to be 

B(\j) = e~ u 'f((x + u' yt)e“')( — u't,\,0). (26) 


Phys. Plasmas, Vol. 5, No. 1, January 1998 


Ellen G. Zweibel 


249 



Equation (26) shows that the magnetic lengthscale shrinks 
exponentially and the held decays at the same rate. The tield- 
lines are also rotated from the y to the .v direction over time, 
because of the shear in the flow. 

Given any incompressible stagnation point flow of the 
type considered here, is there always an initial magnetic held 
which is not amplified by it? This question is easily an- 
swered in the nondegenerate case, in which all the eigenval- 
ues of U are distinct. In this case, the fluid trajectories [which 
are solutions of Eq. (13)] can be written in the form 

x i {t) = a l j(\ 0 )e M i', (27) 

where the fi j are eigenvalues of U, and the a t j are linear 
functions of Xo- It is clear from eqs. ( 11 ) and (27) that the 
elements of the deformation matrix are 



where the partial derivatives are constants. Therefore, D can 
be written as a sum of constant matrices D '\ each one mul- 
tiplied by an exponential function of time. 


D= 2 (29) 

j 

Suppose some i± l has a positive real part; the condition that 
an initial held B 0 not grow with is 

D (,) - B 0 = 0. (30) 

That is, D 0> must be singular, and B 0 must be in the null 
space of D ( ' J . 

In fact, the D <;) are all singular. This follows from the 
condition |D| = 1. Therefore, in a flow in which only one 
eigenvalue has a positive real part, it is always possible to 
find an initial held which is not amplihed by the flow. If two 
eigenvalues have a positive real part, they will not 

necessarily have overlapping null spaces, and so it may turn 
out that all helds are amplihed by the flow. 

Finally, it is worth estimating the amplihcation of the 
held and Lorentz force during the ideal phase, bearing in 
mind that the estimate is based on rough arguments and that 
there are many special cases. Let the eigenvalue with the 
largest positive real part be fi nwx * and let Re(fX milx ) be the 
real part of juL max , Let jm min be the eigenvalue with the most 
negative real part, and Re(fi mi „) be the real part itself. The 
rate of reduction in lengthscale / for the held is dominated by 
fi min , so /~/ 0 exp (Re(fi min )t). The ideal phase ends at the 
time /, when the diffusive timescale becomes equal to the 
flow timescale. 


\/ U ,) tftow : 6 


2 Re(n„ 


-ln (\/ 0 2 t flow ). (31) 


The magnetic held is amplihed at the rate exp (Re(jn max t)). 
According to Eq. (31), the held amplitude at the end of the 
ideal phase is 




In the 2DHS flow given by Eq. (1), fi mux - L fi tm - L so 
the amplihcation factor is proportional to X l/2 . In fact, 


since the real parts of the /x must sum to zero, the exponent 
in eq. (32) must lie between — 1/4 and — I, so the held must 
grow by at least a factor of X 1/4 ( yj i/4 ) during the ideal 
phase. Examples of similar scalings for a variety of dynami- 
cally self-consistent 3D reconnection models are given in 
Reference 13. 

Similar arguments can be made for the growth of the 
Lorentz force F L , which scales as B 2 /!. Thus, at the end ot 
the ideal phase, 

F i F 1 flon) ~ + ^min\ (33) 

The amplihcation rate of the Lorentz force lies between X 1 

and X ~ 5/2 . 

IV. SUMMARY AND DISCUSSION 

In this paper we address the problem of time dependent, 
fast magnetic reconnection in highly conducting fluids in the 
kinematic regime, in which the helds are assumed to be so 
weak that Lorentz forces can be ignored. Reconnection is 
assumed to take place at a stagnation point, and to be essen- 
tially a two stage process. In the hrst step resistivity can be 
ignored, and advection of the held by the flow produces 
structure in the held on very small scales. In the second 
stage, magnetic gradients are large enough that resistivity is 
important, and the held decays. The issue is whether the 
buildup of large Lorentz forces (scaling as an inverse power 
of the resistivity), which would vitiate the kinematic assump- 
tion, can be avoided during the hrst stage. There are two 
main results. 

First, it is not difficult to find examples of flows and 
helds in which reconnection takes place at a rate almost in- 
dependent of the resistivity 77 , and for which, if the kinematic 
approximation is self-consistent initially, it remains so. The 
example discussed in Sec. II is an extension of the 2DHS 
model in which Lorentz forces scale as rj~ m . The crucial 
ingredient is a third component of flow, without which the 
situation would be hopeless. It is argued in Sec. Ill that a 
broad class of incompressible stagnation point flows allow 
fast magnetic reconnection. These results are encouraging for 
some astrophysical problems in which the held and flow ge- 
ometry can be controlled, and the helds are weak, and also 
provide some insight into the nature of almost fast magnetic 
reconnection. 

The second result is that if the initial magnetic held ori- 
entation and coordinate dependence are arbitrary, the field 
will generally be amplihed by the flow. It is the helds which 
are not amplihed which must satisfy special conditions, such 
as B () ^ = 0 for the flow given in Eq. ( 8 ). Since the length- 
scales for the helds shrink as the helds themselves grow, the 
Lorentz forces grow even faster than the helds. Without spe- 
cial symmetries or other restrictions, the forces are amplihed 
by a factor scaling as r) where l=£c/^5/2. This result 
implies that the problem of reconnection of a weak, disor- 
dered magnetic held embedded in a turbulent, highly con- 
ducting fluid cannot be solved self-consistently in the kine- 
matic regime. 

Although amplihcation of the held is generally undesir- 
able in reconnection models for the reasons just described, it 
is exactly the effect sought in models of hydromagnetic dy- 
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namos. The dynamo properties of linear stagnation point 
flows have been investigated by Zel'dovich et ai\ u see also 
Reference 15, by methods which overlap those used in this 
paper. These papers demonstrate the growth of magnetic en- 
ergy with time induced by stagnation point flow, including 
randomly varying flow, but do not estimate the growth of 
Lorentz forces. 

It is worth bearing in mind that these results are obtained 
only in the neighborhoods of stagnation points, not for glo- 
bally bounded flows. This shortcoming is common to other 
reconnection models, which attempt to account for global 
effects through the choice of boundary conditions. 16 

It is often assumed that the effective resistivity of turbu- 
lent fluids is much larger than the Coulomb value, because 
action of the flow on the field produces small scale currents 
which are rapidly dissipated. This is exactly the process stud- 
ied in this paper, although only for local flow models. While 
there can be no doubt that the field can be reconnected at a 
rate that depends only logarithmically on the resistivity, in 
general the fields are amplified and must eventually affect 
the flow. Thus, the results of this paper do not support the 
general concept of a large turbulent resistivity, but instead 
are consistent with arguments made by others based on nu- 
merical computation of global flows, 17 ” 19 or on analytical 
calculations. 20-22 Reference 19, which discusses amplifica- 
tion of the field by a dynamo, shows that Lorentz forces 
reduce the stretching rate of a flow; because of the incom- 
pressibility condition, they must also reduce the shrinkage 
rate. 

In the reconnection models presented here, the field is 
assumed to be initially weak, and the acceptability of the 
models is judged by whether it remains so. These models 
would be on shaky ground in systems such as stellar coronae 
and the bulk of the interstellar medium, which have strong 
fields and require dynamical reconnection theories. The 
models would be better applied to systems such as stellar 
interiors, accretion disks and the early universe, although 
each of these systems has distinctive features arising from 
other physical conditions. 
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